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We model the cosmic medium as the mixture of a generalized Chaplygin gas and a pressureless 
matter component. Within a neo-Newtonian approach we compute the matter power spectrum. 
The 2dFGRS data are used to discriminate between unified models of the dark sector and different 
models, for which there is separate dark matter, in addition to that accounted for by the generalized 
Chaplygin gas. Leaving the corresponding density parameters free, we find that the unified models 
are strongly disfavored. On the other hand, using unified model priors, the observational data 
are also well described, in particular for small and large values of the generalized Chaplygin gas 
parameter a. 



Among the host of models that have been proposed for dark matter and dark energy over the last years, there are 
unified models of the dark sector according to which there is just one dark component that simultaneously plays the role 
of dark matter and dark energy. The most popular proposal along this line is the Chaplygin gas, an exotic fluid with 
negative pressure that scales as the inverse of the energy density [lfl . This phenomenologically introduced equation of 
state can be given a string theory based motivation [2[ . It has also been generalized in different phenomenological ways 
3|. Another example for a unification scenario for the dark sector is a bulk viscous model of the cosmic substratum 
4|. While the Chaplygin gas model (in its traditional and generalized forms) has been very successful in explaining 
the supernovae type la data [5|, there are claims that it does not pass the tests connected with structure formation 
because of predicted but not observed strong oscillations of the matter power spectrum [6J. It should be mentioned, 
however, that oscillations in the Chaplygin gas component do not necessarily imply corresponding oscillations in the 
observed baryonic power spectrum [7|. 

The generalized Chaplygin gas is characterized by the equation of state 
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For A > the pressure p is negative, hence it may induce an accelerated expansion of the universe. The corresponding 
sound speed is positive as long as a > 0. Recently, a gauge- invariant analysis of the baryonic matter power spectrum 
for generalized Chaplygin gas cosmologies was shown to be compatible with the data for parameter values a w 
and a > 3 @. This result seems to strengthen the role of Chaplygin gas type models as competitive candidates for 
the dark sector. Thepresent work provides a further investigation along these lines. While we shall rediscover the 
mentioned results of [g], albeit in a different framework, we also extend the scope of the analysis in the following sense. 
The authors of Q have shown that Chaplygin gas cosmologies are consistent with the data from structure formation 
for certain parameter configurations. Here we ask additionally, whether or not the data really favor generalized 
Chaplygin gases as unified models of the dark sector. By leaving the density parameters of the Chaplygin gas and 
the non-relativistic matter component, respectively, free, we allow for a matter fraction that can be different from the 
pure baryonic part. This is equivalent to hypothetically admit the existence of an additional dark matter component. 
In other words, we do not prescribe the unified model from the start. Moreover, our study is not restricted to the 
spatially flat case. 

Our study relies on a neo-Newtonian approach which represents a major simplification of the problem. In some 
sense, the neo-Newtonian equations can be seen as the introduction of a first order relativistic correction to the usual 
Newtonian equations [l3| . The neo-Newtonian equations for cosmology fl2l. Il3l. [Til fl5j modify the Newtonian equations 
in a way that makes the pressure dynamically relevant already for the homogeneous and isotropic background. This 
allows us to describe an accelerated expansion of the Universe as the consequence of a sufficiently large effective 
negative pressure in a Newtonian framework. While the neo-Newtonian approach reproduces the GR background 
dynamics exactly, differences occur at the perturbative level. However, the GR first-order perturbation dynamics 
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and its neo-Newtonian counterpart coincide exactly in the case of a vanishing sound speed On small scales one 
expects the spatial pressure gradient term to be relevant and the difference to the GR dynamics should be of minor 
importance. Since the observational data correspond to modes that are well inside the Hubble radius, the use of a 
Newtonian type approach seems therefore adequate. 

On this basis our analysis extends previous neo-Newtonian studies to the two-component case. One of the com- 
ponents is a generalized Chaplygin gas, the other one represents pressureless matter. The advantage of employing 
a neo-Newtonian approach is a gain in simplicity and transparency. Our neo-Newtonian approach reproduces the 
parameter estimations for the unified dark matter/dark energy in [8| also numerically. Backed up by this success of 
the neo-Newtonian approach we then enlarge the scope of our analysis and test the validity of the unified model itself 
by relaxing the unified model priors used in Q. Denoting the present value of the Chaplygin gas density parameter 
by r2 c o, we admit the total present matter density parameter r2 m0 to be the sum of an additional dark matter com- 
ponent with density parameter Qdma an d the baryon contribution O;,o, i.e., £l m o = ^dmO + ^fco- Leaving the density 
parameters free, we investigate whether or not the unified model with fl C Q « 0.96, f^o ~ 0.04 and QdmO ~ is favored 
by the large-scale structure data. We mention that a similar investigation using supernova type la data reveals that 
the unification scenario is the most favored one [B[ . 

In the framework of the neo-newtonian formalism, in the conservation equation one takes into account the work 
done by the pressure during the expansion of the universe. At the same time, the equation for the gravitational 
potential must be modified in order to render the equations compatible. This has been done in references (l2l. [T3L [l4| . 
The final equations are 

^ + V-(pv)+pV-v = , (2) 

g + *-w= -5L-v*, (3) 

at p + p 

V 2 = AttG{ P + 3p) . (4) 

For the case of two non-interacting fluids with energy densities p c and p m and pressures p c and p m = 0, respectively, 
the equations are: 

% + V • (p c v c ) + p c V ■ v c = , (5) 
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+ V • (p m U m ) = , (7) 

v m ■ Vu m = -V0 , (8) 



dt 

V 2 <f> = AwG{ Pm +p c + 3 Pc ) . (9) 

The subscript m stands for pressureless matter and the subscript c for the (generalized) Chaplygin gas component. 
Considering now an isotropic and homogeneous universe with p = p(t), p — p(t) and v = ^r, we find 
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Let us define the fractional density contrasts 



S c = ^ and 5 m = 6 -^ (12) 
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for the Chaplygin gas and matter components, respectively, the first-order perturbation equations for the system 
(H-© are 
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where v 2 = j^f- and w c = — . The quantity k 2 denotes the square of the comoving wave vector. Dividing the equations 

([13]) and (| 14[) by Hq and redefining the time as tHo — > t, these equations become dimensionless. In terms of the scale 
factor a as dynamical variable, the system (Ti~3")) - (fr4"|) takes the form 
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where the quantity k 2 denotes the square of the comoving wave vector and Ijj = cH 1 is the present Hubble radius. 
The prime denotes a derivative with respect to a and the definitions 
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have been used. Recall that tt m o = fldmo + ^60- For the unified model to be an adequate description one expects 
il m o ~ tlbo- In case the data indicate a substantial fraction of fldmo, the unified model will be disfavored. 
The power spectrum is defined by 



■P = si 



(22) 



where 8k is the Fourier transform of dimensionless density contrast 8 m . We will constrain the free parameters using 
the quantity 
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where V° is the observational value for the power spectrum, V\ is the corresponding theoretical result and Oi denotes 
the error bar. The index i refers to a measurement corresponding to given wavenumber. The quantity (|23[) qualifies 
the fitting of the observational data for a given theoretical model with specific values of the free parameters. Hence, 
X 2 is a function of the free parameters of the model. The probability distribution function is then defined as 



F(x n )=F e-x 2 ^/2 



(24) 



where the x n denote the ensemble of free parameters and Fq is a normalization constant. In order to obtain an 
estimation for a given parameter one has to integrate (marginalize) over all the other ones. For a more detailed 
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FIG. 1: The two-dimensional probability distribution function (PDF) for fldmO an d ^AO (left) and the corresponding one-dimensional 
probability distribution functions for the non-flat ACDM model. In the left panel: the darker the color, the smaller the probability. 




FIG. 2: The results for the general case with four free parameters. From left to right: the one-dimensional PDFs for a, A, Q c o and f^dmO- 

description of this statistical analysis see reference [|| . From now on we focus on the 2dFGRS observational data for 
the power spectrum [TB|. We use the data that are related with the linear approximation, that is, those for which 
k h^ 1 < 0.185-Mpc -1 , where h is defined by Hq = 100 • ftkm/s • Mpc. This definition should not be confused with the 
preceding definition of the function h(a). To fix the initial conditions we use the BBKS transfer function [9(. This 
procedure is described in more detail in references [lol [Tl| . 

To "gauge" our approach, let us first consider the ACDM model. In the general (non-flat) case there are two 
parameters: Qdmo an d ^ao- hi figure [T] we show the two-dimensional probability distribution function (PDF) as well 
the one-dimensional PDFs for the dark matter parameter Sldmo and for the cosmological constant parameter Qao, 
respectively. From the two dimensional graphic it is clear that there is a large degeneracy for the parameter Qao, 
while the region of allowable values for Qdmo lS quite narrow. The degeneracy for the cosmological constant density is 
less visible in the one-dimensional PDF graphic, but it is still considerable. Incidentally, the minimum value for the 
X 2 parameter is 0.3822 for fld m a = 0.2387 and f^Ao = 0.5937, corresponding to an open universe. 

The four free parameters to be constrained in our Chaplygin gas model are ftdmo, ^co, A and a. The one-dimensional 
PDFs for a, A, fldmo and O c o are displayed in figure [H It can be seen that the preferred values are either a C 1 or 
a > 2, while the probability is higher for large values of fldmo and small values of £l c o. This show clearly that the 
unification scenario is disfavored. 



If the unification scenario with dark matter and dark energy as a single fluid in a spatially flat universe is imposed 
from the beginning, the results of reference [|[ are essentially confirmed: there are parameter ranges for which the 
data are well described by the generalized Chaplygin gas model, see figure |3] The probability distribution function 
for a is high for very small (near zero) or very large (greater than 2) values of a. Allowing the parameter A to vary, 
we find that its one-dimensional PDF initially decreases with A, but increases as A = 1 is approached. Notice that 
values a > 1 imply a superluminal sound speed and are therefore unphysical (see, however, [8[). 

What is the origin of these apparently contradictory results? The first aspect to be mentioned is that the matter 
power spectrum data only poorly constrain the dark energy component. Even for the ACDM model the matter power 
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FIG. 3: The results for the flat case with Q b0 = 0.043, f^mO = and Q c q = 0.957, corresponding to the unification scenario. From left 
to right: the two-dimensional PDF for a and A (with the same color convention as before), the best fitting curve for the power spectrum, 
the one-dimensional PDFs for a and A. 



spectrum gives information mainly on the dark matter component, the dark energy component remaining largely 
imprecise. It is not by chance that the dark energy concept emerged from the supernova data. Our results for the 
Chaplygin gas model show that a large amount of dark matter, different from those described by the Chaplygin gas, is 
necessary to fit the data. However, the dispersion is quite high. For the flat case with a three-dimensional parameter 
space we find at 2er, that fi^mo = l^o gi- Another point is the use of the neo-Newtonian formalism. However, for 
small values of the parameter a, the main case of interest here, the differences to the full general relativistic treatment 
are not expected to be substantial. Moreover, in the cases of overlap the results of the full theory are reproduced. 
Finally, possible statistical subtleties may influence the outcome of the investigation. But as far as we could test the 
statistical analysis (precision, crossing different information, etc), the results seem to be robust. If this is really the 
case, we must perhaps live with the fact that, while the SNe type la data favor a unified model of the dark sector [f|, 
this scenario is disfavored if large scale structure data are taken into account, unless specific priors are imposed. 
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